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Superfluidity in fermionic systems originates from pairing of fermions, and Bose condensation of 
these so-called Cooper pairs. The Cooper pairs are usually made of fermions of different species; 
thus the most favorable situation for pairing and superfluidity is when the two species of fermions 
that form pairs have the same density. This paper studies the possible superfluid states when the 
two pairing species have different densities, and show that the resultant states have remarkable 
similarities to the phases of liquid crystals. This enables us to provide a unified description of the 
possible pairing phases, and understand the phase transitions among them. 



In the Bardeen-Cooper-Schrieffer (BCS) theory for 
fermionic superfluidity, fermions form Cooper pairs that 
condense into the zero momentum state. A Cooper pair 
is often made of fermions of different species; in super- 
conductors they are electrons of opposite spins. Thus 
the most favorable situation for pairing is when the 
two species of fermions have the same density, so that 
there is no unpaired fermion in the ground state. It 
has been a long-standing fundamental question as to 
what kind of pairing states fermions can form when the 
two fermion species have different densities or chemi- 
cal potentials. An early suggestion was due to Fulde 
and FerrellQ], and Larkin and OvchinnikovQ, who ar- 
gued that the Cooper pairs may condense into a finite 
momentum state; this Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state breaks translation and rotation symme- 
tries. More recently other suggestions have been put for- 
ward, including deformed Fermi surface pairing^, 0, 
and breached pairing[(J 0, EJ states, each with their dis- 
tinct symmetry properties. Here we show that all these 
states arc different phases of a quantum liquid crystal. 
This allows us to provide a unified description of all these 
phases, propose a global phase diagram, and understand 
the nature of the transitions between neighboring phases. 

The issue of pairing between unbalanced fermion 
species originally arose in superconductors subject to ei- 
ther an external magnetic field or internal exchange field, 
which gives rise to Zeeman splitting between up- and 
down-spin electrons that form Cooper pairs. More re- 
cently the same issue has been under close scrutiny in the 
context of pairing and superfluidity in nuclear matter 4 , 
neutron stars0, and high density quark matter[^,0, !j, 
where the unbalance is due to more intrinsic effects like 
difference in mass of the fermions that form pairs. Per- 
haps the most promising place where some of these novel 
phases can be directly observed are trapped cold atom 
systems jl2j , where one has the best control of the 

properties of the constituents and strength of pairing in- 
teractions. Thus this is a fundamental issue that is of 
importance to all branches of physics. To better eluci- 
date the properties of and relations between the various 



proposed phases, in this work we will use concepts devel- 
oped in studies of liquid crystals, which is traditionally 
viewed as a branch of classical physics. The analogy to 
liquid crystals allows us to determine the global phase 
diagram of the system, and gain insight into the nature 
of the transitions between neighboring phases. 

We start our discussion with the FFLO state, which 
has the longest history of studies, and very strong exper- 
imental evidence for its existence has been found recently 
in a heavy fermion superconductor, CeCoIn^ [l3l [p3 | . Fol- 
lowing the superconductivity terminology, throughout 
this paper we will use "spin" indices a =1\ I to label the 
two different species of fermions that form Cooper pairs, 
"Zeeman splitting" A/i = /it — fi^ to represent the chem- 
ical potential difference between the two species, and 
"magnetization" m to represent their density difference. 
When A/i ^ 0, up- and down-spin electrons form Fermi 
seas with different Fermi momenta _p ft and pfi in the 
normal state; it was thus suggested [ll|2| that when pair- 
ing interaction is turned on, the electrons with opposite 
spins on their respective Fermi surfaces would pair up to 
form a Cooper pair with a net momentum p « pp-\ —Pfi- 
This results in a pairing order parameter A(r) that is os- 
cillatory in real space, with period 2-Kh/p. In general the 
structure of A(r) is characterized not just by a single 
momentum p, but also by its higher harmonic compo- 
nents. More detailed mean-field study 15] suggested the 
following real space picture for FFLO state: it is a state 
with a finite density of uniformly spaced domain walls; 
across each domain wall the order parameter A (which 
is real in the mean-field theory) changes sign, and the 
excess magnetization due to spin unbalance are localized 
along the domain walls, where A (which is also the gap 
for unpaired fermions) vanishes; see Fig. la. Thus the 
total magnetization is proportional to the domain wall 
density. This picture was made more precise by an exact 
solution in one-dimension (ID) based on bosonized de- 
scription of spin-gapped Luttinger liquids ^(J, where the 
domain walls arc solitons of the sine-Gorden model that 
describes the spin sector; each soliton carries one half- 
spin. While quantum and thermal fluctuations do not 
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FIG. 1: (color online). Illustrations of various possible pair- 
ing states, a, Domain walls in the mean-field description 
of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, where 
the pairing order parameter (or gap) changes sign and the 
magnetization is located, b, Domain walls of the FFLO state 
are replaced by chains of alternating vortices and antivortices 
once fluctuations are taken into account. The red and blue 
lines are where the real and imaginary parts of the order pa- 
rameter vanish respectively, and their intersections are the 
vortices. The ordered chains break translation and rotation 
symmetry in a manner identical to the smectic phase of a liq- 
uid crystal, c, The nematic phase of the vortex/antivortex 
chains (or their segments, represented by black lines), which 
correspond to the deformed Fermi surface pairing (DFSP) 
state, d, The isotropic phase of the vortex/antivortex chains 
(or their segments), which correspond to the breached pairing 
(BP) state, e, The BCS state in which vortices and antivor- 
tices (generated by thermal or quantum fluctuations) form 
closely bound pairs, f, The normal state in which the vortices 
and antivortices are unbound. 



allow true long-range order in ID, such order can be sta- 
bilized by weak interchain couplings[16| . Coming back to 
isotropic high D cases, it is clear that the presence and 
ordering of these domain walls break rotation symmetry, 
and translation symmetry in the direction perpendicular 
to the walls, although translation symmetry along the 
wall remains intact. Thus the symmetry properties of the 
FFLO state is identical to that of the smectic phase of 
liquid crystals (smectic- A phase to be more precise) 01 • 
The domain wall picture described above needs to be 
modified when fluctuations above mean-field theory are 
taken into account. In mean-field theory the pairing or- 



der parameter A is real; it thus vanishes along domain 
walls which are lines in 2D and sheets in 3D. In the pres- 
ence of fluctuations however A is complex, thus it only 
vanishes where the real and imaginary parts of A are 
simultaneously zero (see Fig. lb); these points in 2D 
and lines in 3D are nothing but the familiar vortices. 
Thus in the presence of fluctuations the mean-field do- 
main walls are unstable and replaced by chains that are 
made of alternating vortices and antivortices (see Fig. 
lb), and the FFLO state is a state in which these chains 
are generated by the Zeeman splitting, and they line-up 
to break the translation symmetry perpendicular to the 
chain direction[l8|. 

The appearance of these vortex-antivortex chains when 
the fermionic superfluids are subject to a sufficiently large 
Zeeman splitting can also be understood from the fol- 
lowing consideration. Zeeman splitting tends to gen- 
erate spin-polarized quasiparticles; in a uniform super- 
fluid state this can happen only when the splitting is 
large enough to overcome the quasiparticle gap. On the 
other hand it is known [l(il| that the quasiparticle excita- 
tions are gapless inside the vortex core; thus the vor- 
tices are easily polarizable, making their presence en- 
ergetically favorable when Zeeman splitting is present. 
The net vorticity has to remain zero so that equal num- 
bers of vortices and antivortices must be generated, and 
since vortices repel/attract each other when they have 
the same/opposite vorticity, they naturally form an alter- 
nating pattern. This naturally leads to the vortex chain 
configurations discussed above. 

When quantum and/or thermal fluctuations are suf- 
ficiently weak, these infinitely-long vortex/antivortex 
chains line up and form an ordered smectic phase that 
breaks translation symmetry along the direction perpen- 
dicular to the chains and rotation symmetry, as illus- 
trated in Fig. lb; this corresponds to the FFLO state in 
the presence of fluctuations. In classical liquid crystals 
it is known that as one increases thermal fluctuations, 
the broken symmetries of the smectic phase are restored 
in the following sequence ^3 : the translation symmetry 
is restored first when the smectic melts into a nematic 
that breaks the rotation symmetry only, and then the 
nematic melts into an isotropic liquid that has no bro- 
ken symmetry. We thus expect the same sequence of 
phases and phase transitions occur in superfluids with 
unbalanced fermion pairing, as we increase the strength 
of either thermal or quantum fluctuations. More specif- 
ically, as fluctuations increase, these infinitely long vor- 
tex/antivortex chains can break into segments so that 
dislocations can appear in the lattice formed by these 
chains; when such dislocations proliferate and unbind the 
translation symmetry is restored, although (segments of) 
the chains are still aligned in a preferred direction, so that 
the rotation symmetry remains broken (see Fig. lc for 
an illustration). This is the analog of the nematic phase, 
which has the same symmetry as the deformed Fermi sur- 
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face pairing (DFSP) stateHlill; we thus identify the 
DFSP state as the nematic. The chains (or their seg- 
ments) will lose their orientation alignment upon further 
increasing fluctuations (see Fig. Id for an illustration), 
thus restoring the rotation symmetry (through prolifer- 
ation of disclinations); this isotropic liquid phase is thus 
identified with the breached pairing (BP) phaseHHil. 
While all broken spatial symmetries are restored by fluc- 
tuations at this point, these phases are still superfluid 
phases with a spontaneously broken U(l) gauge symme- 
try, as long as the vortices and antivortices are bound 
together to form chains or their segments. Further in- 
creasing fluctuations will eventually break the segments 
apart into unbound vortices; this is the normal phase in 
which the U(l) gauge symmetry is also restored (see Fig. 
If). 

It is clear that the fundamental similarity among the 
three different superfluid phases discussed above is the 
appearance of the vortex/antivortex chains (or their seg- 
ments) in the ground state; they are responsible for ac- 
commodating the excess magnetization or unbalance be- 
tween the fermion species. This is very different from the 
BCS phase. While in the BCS phase (or ordinary super- 
fluids due to simple Bose condensation as described by 
the XY model) thermal or quantum fluctuations can still 
generate vortex excitations, they appear as bound pairs 
each made of a vortex and an antivortex (see Fig. le), 
instead of forming chains or segments of chains. 

Based on the discussion above, we propose a phase di- 
agram in 3D (Fig. 2) for pairing between unbalanced 
fermion species, assuming the Zeeman splitting A/x is 
strong enough to destabilize the BCS state, yet not strong 
enough to eliminate pairing and superfluidity at T = 0. 
For pairing fermions with the same mass, this requires 
Afi ~ Aq, where Ao is the pairing gap at T = and 
Afi = 0; it parametrizes how much energy it costs to 
break a Cooper pair. We use Ao as our unit of en- 
ergy (although it depends on the pairing strength and 
other parameters itself). As usual thermal fluctuation 
is controlled by temperature. Quantum fluctuation is 
controlled by the strength of pairing interaction; it can 
be measured by the dimensionless ratio A /Ep, where 
Ep = pp/2M is the Fermi energy in the absence of A/i 
and pairing interaction, and M is the fermion mass. In 
trapped cold fermionic atom systems the pairing strength 
can be controlled by tuning the distance to a Fesh- 
bach resonance by a magnetic field, where one can ex- 
plore a wide parameter range from the weak coupling 
(or BCS) to strong coupling (which corresponds to the 
Bose-Einstein condensation or BEC of diatom molecules) 
regimes of superfluidity [2(J. As we increase thermal (in- 
creasing T) and/or quantum (increasing Ao/ Ep) fluc- 
tuations, the smectic phase (that represents the FFLO 
state) first melts into a nematic (representing the DFSP 
state), and then yields to the isotropic (representing the 
BP state). The experimental detection of these phases 
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FIG. 2: (color online). Schematic phase diagram for pairing 
between unbalanced fermion species in 3D. The FFLO state, 
which is a smectic, is stable at low temperature and weak cou- 
pling. As thermal fluctuation (controlled by T) and quantum 
fluctuation (controlled by pairing strength and parametrized 
by the ratio between pairing gap and Fermi energy, Aq/ Ef) 
increase, it first melts into a nematic (or the deformed Fermi 
sea pairing state), which then melts into an isotropic liquid 
crystal (or breached pairing state). Direct transitions between 
these superfluid states and the normal state are also possible 
at high temperature. 



should be fairly straightforward in these svstems[nlll2j]. 
The nematic-isotropic transition is first order due to a 
cubic term in the appropriate Landau theory |l7j. The 
smectic (-A)-nematic transition is (typically) a continuous 
transition, although the experimentally observed critical 
behavior is not yet completely understood^]- Direct 
transition between these superfluid states and the nor- 
mal state are also possible. The transition between the 
BP and normal phases is in the same universality class 
as the usual superfluid (or XY) transition, which is usu- 
ally continuous. On the other hand the FFLO/DFSP 
to normal transitions should be first order even if the 
mean-field theory suggests a continuous transition; this 
is because both spatial and gauge symmetries are broken 
simultaneously here. Such fluctuation-driven first-order 
transition has been explicitly demonstrated for the FFLO 
case recently j^. The situation is quite different in 2D; 
there the smectic phase is unstable for any finite temper- 
ature, while the nematic-isotropic transition is expected 
to be of the Kosterlitz-Thouless type 22]. On the other 
hand the quantum (or T = 0) phase transitions among 
these phases are poorly understood at this point, and will 
be the subject of future work. 

In the following we present a concrete case of the 
breached pairing (BP) state, in the strong-coupling (or 
BEC) regime. In this case unpaired, spin-polarized 
fermions co-exist with closely-bound fermion pairs (or 
molecules) that Bose condense; the appropriate fermionic 
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many-body wave function takes the form: 



i*) = n it n <«* 



a)l°>> (!) 



iki<fe' 



ik|>fe' 



where |0) is the vacuum state, k' F 



[6vr 2 (n T - n x )]V3 
is the Fermi wave-vector of the excess (or unpaired) 
fermions, and u k and v k (which satisfy u 2 . + = 1) 
are coherence factors; they are determined by identifying 
<h k = V]c I uy as the (momentum space) molecular wave 
function|23j|. which in turn can be obtained by solving 
the two-body Schrodinger's equation with the constraint 
that $ k = for |k| 



<k' F : 



h 2 k 7 
M 



*k + ^^kk'*k' = -2|A|$ k , 



(2) 



where Vkk' is the pairing potential in momentum space. 
This only leads to a very minor modification of $ k from 
the genuine two-body problem in the BEC limit, be- 
cause components with |k| < k' F only carry a very small 
weight for the closely bound state with size much smaller 
than l/k' F ; thus the pairing is only "breached slightly" 
by the presence of the spin-polarized, unpaired fermions. 
Clearly this is the energetically favored state at strong 
coupling (compared to the FFLO and DFSP states), as 
it minimizes the kinetic energy of the unpaired fermions 
while preserve the pairing correlation for the remaining 
fermions that do form pairs. The wave function takes 
a form similar to the one proposed by Liu and WilczekQ; 
the main difference is that here the unpaired fermions oc- 
cupy states near the origin in momentum space (instead 
of near fc^ of the non-interacting Fermi gas 6J) and form 
a circular Fermi sea of itself; the origin of the difference is 
we are considering the strong coupling case while Liu and 
Wilczek were considering the weak-coupling limit. But 
they share the same symmetry property (both isotropic), 
and are therefore different forms of the BP phase. 

We close by noting that quantum liquid crystal 
phases have also been proposed to describe other 
strongly interacting fermionic systems, including; cuprate 
superconductors |2J| and quantum Hall liauids |25ll26| . 
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This work was supported by by National Science Foun- 
dation grant No. DMR-0225698. 

Note Added - After the present paper has been sub- 
mitted for publication, the author became aware of three 
interesting new Dreprints |27l 1281 129| that address closely 
related issues from somewhat different angles. 
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